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APPENDIX  A  - 1988  DOPPLER  SURVEY  OF 
AFIT BUILDING  640  ROOFTOP 


DEFENSE  MAPPING  AGENCY 

hydrographic/topographic  center 

TERRESTRIAL  COMPUTATIONS  BRANCH 
UASHIN6T0N,  DC  20315-0030 


PROJECT:  ANTENNA  SURVEY 


LOCATION:  MRIGHT-PATTERSON  AFB,  OHIO 


DATUMS: 


GEOID  HEIGHT 
(H«h-H). 


1.  WORLD  GEODETIC  SYSTEM  OF  1972  (HGS  72) 

2.  WORLD  GEODETIC  SYSTEM  DF  1984  (WGS  84) 

3.  ELEVATION  BASED  ON  WGS  84  DERIVED  ELLIPSOID  HEIGHT  AND 
COMPUTED  6RAV1METRICALLY  USING  MGS  84  POTENTIAL  COEFFICIENTS 


SURVEY  DATE:  JULY  1988. 

1.  GEODETIC  POSITIONS,  ELLIPSOID  HEIGHTS  AND  DOPPLER  DERIVED  ELEVATIONS. 


NORTH 

STATION  DATUM  LATITUDE 
DEG  MIN  SEC 


DOPP  32058 

WGS 

72 

39 

46 

55.306 

WGS 

84 

39 

46 

55.424 

ALGN  MARK 

WGS 

84 

39 

46 

55.128 

GPS  ANT. 
(BOLT) 

WGS 

84 

39 

46 

55.369 

GPS  ANT. 

(TOP  OF  BASE 
OF  MOUNT) 

LOADING  DOCK 


WEST 

LONGITUDE 

DEG  MIN  SEC 

ELLIPSOID 

HEIGHT 

METERS 

ELEVATION 

METERS 

084  04  58.472 
084  04  57.918 

269.41 

271.77 

306.27 

084  04  59.648 

272.42 

306.92 

084  04  59.422 

272.54 

307.04 

311.84 

293.67 

•ALTHOUGH  ELEVATIONS  ARE  GIVEN  TO  0.01  METERS,  Jlnm 

AND  HAVE  AN  ACCURACY  OF  ♦/-  5  METERS  DUE  TO  THE  UNCERTAINTIES  OF  THE  GEOID 

HEIGHT. 
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SUMMARY  OF  S  A  T E L L I T E - 0 B S E R V E D  STATION 


UOCATIOM  OOFFVtR  NO. 

WRIGHT-PATTERSON 

OAYTOM,  OHIO  32058  • 

ITAMRINO  ON  MANN 

OOPP  STA  32058  1988  GSGS  ' 

ACCHCT  ICAST  IN  HARKt 

TTVC  or  STATION  KAMA 

OMA 

BRONZE  DISK 

tOUINMCNT/fCNIAL  HO. 

MAG  1502-05*688 

HEIGHT  OF  TNACKIMC  EOOIPMCMT  NCFCNCNCE 
NOIMT  ABOVE  STATION  mark:  2.412m  ’ 

TNACKINO  EQUIPMENT  REFERENCE  POINT 

RED  BAND  ON  ANTENNA 

OBSCNVeO  BY  lACCNCYl 

OHAHTC/GSGS 

satclutcisi  OOSCNVCO  30130, 

30700.  30240.  30480.  30500  ’ 

PCNIOO  or  OCCUPATION 

07  -  10  JULY  1988  ’ 

iATgLUTE-DERIVgO  COO^DIMATES 


passes  Acecpfco 

pcenccs  of  • 

MESIOUAL  RMS 

STATION  SET 

ONAVITV  model 

ELLIPSOID 

MINIMUM  CLEV. 
ANCLE:  ^ 

48 

mecoOM: 

0.09m 

WGS  84 

WGS  84 

WGS  84 

5  OEG. 

(Softtii9<‘d*n^tdco9rdin9itt  r^Urrmd  to  gtoHon  mark) 


H  39  46  55.424 


SQ6  020.04m 


E  275  SS  02.082 


-4  6B2  282.48m 


271.77m 


4  059  S93.06m 


AccuwAcr  • 

1.5  METERS  IN 
EACH  AXIS  (90X 
LINEAR  ERROR) 


ccofdlnofwt  of  $foH0ft  mark  teontfemd  <o 


WGS  72) 

DATUM 


H  39  46  S5.306_ 


E  27S  SS  01.528_ 


269.41m 


WGS  72 
ewuiMoio 


506  006. e2m 


-4  882  282.76m 


4  0S9  587. 68m 


WGS  72 

DATE  OF  TMANSFOAMATION 


AY 


AZ 


A 

DATUM  (HORIZONTAL) 

ELLIPSOID 

iATC  OF  AOajSTMSnT  ORDER 

SURVEY  BY  lACEHCYl  DATE  LOCATION  OF  1 

lURVEY  DATA 

elevation  imi  H^h-N 

306.3m 

DATUM  (vertical) 

GEOIO  HEIGHT  (N)  (WGS  84) 
-34..5m  t  Smdc) 

ELLIPSOID  HEIGHT  (h) 

ONOCP  (CLEv.t  established  BY  1 

AGENCY)  DATE 

SOUHCCOTINI  COMPUTED  GRAV I  METRICALLY 

USING  WGS  84  POTENTIAL  COEFFICIENTS 

FROM 

TO  1 

(  )  AS  FROM  north 

DISTANCE 

DOHP  32058 

267  18  01.3 

35.839m 

UUPP  32Ub8““ 

aLgN  HAkk 

Zb7  30  53.4 

42.159m 

Hit  ST  ATION 

STATION  OCCUPATION  MCPOMT 


CCOCCTtC  iNrOAMATlOH  WEPQWT 


STATION  DCSCAIPTION 


Sunvcy  DIACMAM 


STATION  SITE  SKETCH 


PNOTOIOEMnnCATIOM 


ASTRONOMIC  COOAOl NATES 


station  photos 


LOCATIQH/WEMARKS 

DMAHTC _ 


UMAHTC 


umrnr 


UHAHTC 


AEPARCO  SY/OATC 

OMAHTC/RBX/OCT  88 


CHECKED  BV/OATC 

OMAHTC/BHS/OCT  88 


revised  by/oate 


CHECKED  BY/OATC 


DMA  FORM  IZtO.I.R 
*C»  >• 
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fOI 


CPS-v%< 

AuTCmM^ 


/MoBlt-E 

necei'/tt 

AiJTEMMA 

BAse 

/^  S-rflTiotJ  ,. 

AtAM 

MAAK 


SToBACC 

ILOOA\ 


LEGEND 

▲  .  TRACKING  STATION 
O  -  SURVEY  STATION 


-  PHOTO  POINT  A  -  GEODETIC  STATION/RM'S  t  BN* 

•  I  ^  .f' 

-  MEASURED  DISTANCE  !  / 


APPENDIX  B -LATITUDE/LONGITUDE/ALTITUDE 
TO  ECEF  TRANSFORMATIONS 


The  Earth  Center  Earth  Fixed  (ECEF)  coordinate  system  is  commonly  used  in 
GPS  calculations.  In  the  ECEF  system,  the  positive  'x'  axis  is  oriented  in  the  direction 
of  zero  degrees  longitude,  the  positive  'y'  axis  points  in  the  direction  of  90  degrees 
East  longitude.  The  x-y  plane  is  coincident  with  the  Earth's  equatorial  plane  and 
rotates  with  the  earth.  The  positive  'z'  axis  is  normal  to  the  x-y  plane.  For  a  detailed 
explanation  of  reference  coordinate  systems  the  reader  is  referred  to  [KAP96]. 

Satellite  ephemeris  data  is  provided  in  the  Earth-Centered  Inertial  (ECI) 
coordinate  system.  These  parameters  are  converted  through  transformation  matrices 
to  the  ECEF  coordinate  system.  Transformation  from  the  ECEF  system  to  latitude, 
longitude  and  altitude  normally  completes  the  process  for  the  end  user.  For  the 
purpose  of  this  thesis,  the  ECEF  position  of  the  MobUe  XR5-M  receiver  was 
compared  to  the  ECEF  position  of  the  XR5-M  Base  Station. 

One  of  the  first  challenges  was  to  determine  a  reasonably  accurate  position  for 
the  Base  Station.  Since  the  primary  measurement  of  interest  was  the  relative 
difference  between  the  Base  Station  antenna  and  the  Mobile  antenna,  the  Base  Station 
was  used  in  a  relative/differential  GPS  configuration.  By  using  relative/differential 
GPS,  most  of  the  position  error  attributed  to  Selective  Availability  (SA), 
ionosphere/troposphere,  satellite  perturbations,  ephemeris  prediction  and  satellite 
clock  bias  were  removed.  The  major  sources  of  position  error  between  the  Base 
Station  and  the  Mobile  receiver  was  a  result  of  multipath,  receiver  measurement  noise 
and  RF  interference.  This  setup  was  used  to  allow  for  a  reasonable  method  to  measure 
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the  effects  of  jamming.  The  Mobile  receiver  antenna  was  physically  located  4.7 
meters  approximately  east  of  the  Base  Station  antenna.  Position  ECEF  coordinates  of 
the  Mobile  receiver  antenna  were  corrected  in  MATLAB®  software  during  post¬ 
processing  back  to  the  position  of  the  Base  Station  antenna  so  that  the  antennas  would 
be  at  the  "same"  position.  Position  error  was  determined  by  taking  the  difference 
between  the  corrected  Mobile  antenna  location  as  measured  by  the  Mobile  receiver 
and  the  "ALGN  MARK"  position  of  the  previous  1988  survey  (Appendix  A).  No 
adjustment  was  made  to  correct  the  actual  coordinates  of  the  Base  Station  to  the 
surveyed  location  of  the  "ALGN  MARK"  because  the  relative  position  between  the 

two  antenna  was  the  measurement  of  interest. 

The  rooftop  of  Building  640  was  surveyed  by  the  Defense  Mapping  Agency  in 
1988;  details  of  which  are  provided  in  Appendix  A.  Doppler  Station  32058  was  the 
main  survey  location  and  two  additional  surveyed  sights  were  provided,  ALGN 
MARK”  and  “GPS  ANTENNA  (BOLT).”  For  the  minimum  cost  and  ease  of  antenna 
installation,  the  Base  Station  antenna  was  located  approximately  8.2  meters  east  of  the 
previous  surveyed  “ALGN  MARK"  position  fi-om  the  1988  survey.  The  survey 
reported  an  accuracy  of  1.5  meters  in  each  axis  (90  %  linear  error).  Elevation  was 
given  as  306.92  meters  but  because  the  heights  were  Doppler  derived  the  elevation 
measurement  had  an  accuracy  of  ±5  meters. 

The  values  for  the  "ALGN  MARK"  position  were  used  as  coordinates  for  the 
Base  Station.  The  geodetic  coordinates  for  “ALGN  MARK”  taken  fi'om  the  1988 
survey  were  determined  to  have  the  following  properties: 
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Lat  =  39  46  55.128  N 
Long=  084  04  59.648  W 

Ellipsoid  Height  =  272.42  meters,  Altitude  306.92  meters. 

On  the  XR5-M  receiver,  input  of  the  above  parameters  on  the  computer  display 


initialization  page  resulted  in  the  following: 


Lat  39  46  55. 128  N  was  displayed  as  39  46.9188  N 
Long  084  04  59.648  W  was  displayed  as  084  04.9941  W 
Initial  Altitude  :  306.92  was  displayed  as  306.9 


The  XR5-M  variable  “npposn”  provides  the  receiver's  best  estimate  of  position 
in  ECEF  coordinates.  When  the  unit  is  configured  as  a  Base  Station  the  value  input 
for  latitude,  longitude,  and  altitude  is  fixed  to  the  values  used  in  the  initialization  menu 
for  the  Base  Station  location.  The  coordinates  in  ECEF  (meters)  for  "ALGN  MARK" 
as  converted  by  the  XR5-M  firom  the  above  input  were  (in  meters). 

X=  505979.8731834281 

Y  = -4882293.899751916 

Z=  4059587.157385162 

The  1988  survey  did  not  provide  ECEF  coordinates  for  the  "ALGN  MARK 
position.  As  a  check  on  the  XR5-M  internal  calculations,  the  coordinates  for  the 
"  ALGN  MARK"  position  fi-om  the  1988  survey  were  converted  to  ECEF  using 
MATLAB®  code  and  a  more  accurate  set  of  transformation  equations  fi'om  [KAP96] 
(as  compared  to  those  presented  in  AHT  course  EENG  533).  The  resulting  ECEF 
coordinates  (meters)  for  the  "ALGN  MARK"  were  as  follows: 

X=  505979.7386127784 

Y  =  -4882293.041839992 

Z=  4059586.447032988 
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The  differences  between  the  surveyed  "ALGN  MARK”  ECEF  coordinates  as 
calculated  by  the  XR5-M  receiver  compared  to  using  the  equations  provided  by 
Kaplan  using  MATLAB®  to  calculate  ECEF  position  (in  meters)  were  as  follows. 

X  = -0.1345706497086212 
Y  =  0.8579119239002466 
Z  =  -0.7103521740064025 

A  similar  comparison  was  conducted  using  the  transformation  to  ECEF 
coordinates  provided  in  the  1988  survey  for  the  "Doppler  Station  32058."  Values  for 
ECEF  coordinates  for  the  Doppler  Station  provided  in  the  survey,  the  XR5-M 
calculated  ECEF  coordinates,  and  the  MATLAB®  calculated  ECEF  coordinates  using 
Kaplan's  equations  are  summarized  in  Table  B.l. 


Latitude 
Longitade 
Ellipsoid  Height 
for  Doppler  Station 

ECEF  Coordinates 
from  1988  Survey 
(x,y^) 

in  meters 

ECEF 

transformation  by 
XR5  (x,y,z) 
in  meters 

ECEF 

transformation  using 
Matlab  and  Kaplan's 
equations  (x,y,z)  in 
meters 

N  39  46  55.424 

506020.04 

506020.1624427695 

506020.0340076035 

E  275  55  0X082 

^882282.48 

-4882283.308111367 

271.77^ 

4059593.06 

4059593.728202703 

4059593.047042757 

Table  B.l  -  Conversion  of  Latitude/Longitude/Altitude  to  ECEF  Coordinates 


Differences  between  the  surveyed  “Doppler  Station  32058”  ECEF  position,  the 
XR5-M  calculated  ECEF  position  and  the  MATLAB®/Kaplan  calculated  ECEF 
position  (in  meters)  is  provided  in  Table  B.2.  Based  on  the  calculations  in  Table  B.2 
for  ECEF  position  differences,  the  method  of  using  MATLAB®  and  Kaplan's 
equations  were  more  accurate  than  the  XR5-M  receiver  calculations  (assuming  that  the 
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Defense  Mapping  Agency  calculations  for  ECEF  coordinates  were  still  accurate). 
Based  on  the  results  in  Table  B.2  the  conversion  process  within  the  XR3-M  receiver 
contributed  to  3D  position  error  on  the  order  of  1.07  meters.  The  conversion  process 
using  MATLAB®  code  and  Kaplan's  equations  contributed  to  3D  position  error  on  the 
order  of  approximately  0.0177  meters. 


Difference  between  ECEF 
Coordinates  from  1988 
Survey  and  XR5-M 
transformation  (x,y,z) 
in  meters 

Difference  between  ECEF 
Coordinates  from  1988  Survey 
and  Kaplan/MATLAB^ 
transformation  (x,y^)  in  meters 

-1.224427695269696e-01 

-5.992396443616599e-03 

8.281113663688302e-6n 

-1.042952854186296e-02 

-6.682027028873563e-01 

-1.295724278315902e-02 

Table  B.2  Differences  in  Conversion  to  ECEF  coordinates  for  Doppler  Station 


The  position  calculated  by  the  XR5-M  for  the  "ALGN  MARK"  was  considered 
acceptable  and  was  used  with  the  knowledge  that  it  could  introduce  calculation  errors 
of  one  meter  magnitude  into  final  results  data. 

Any  fiiture  research  work  using  the  Base  Station  receiver  antenna  location  to 
measure  absolute  position  should  be  aware  of  the  8.2  meter  difference  between  the 
“ALGN  MARK”  coordinate  fi’om  the  1988  survey  and  the  actual  Base  Station  antenna 


location. 


APPENDIX  C-XR5-M  SYMBOL  VARIABLES 


A  symbol  file  was  generated  using  the  Data  Monitor  program  provided  with 
the  XR5-M  software.  Prior  to  using  the  symbol  file  <kdsymbl3.sy>  it  was  necessary 
to  load  the  main  symbol  file  name  <ml237.sy>.  For  further  details  the  reader  should 
refer  to  the  XR3-M  Control  and  Display  Unit  Software  Manual,  the  XR5-M  Power 
Users  Data  Monitor  Toolkit  Manual,  and  the  data  dictionary  text  file  <datadic2.txt  > 
provided  with  the  XR5-M  receiver  software  [NAV93,  NAV96].  The  variables  in  this 
Appendix  were  recorded  during  Test  12;  previous  tests  consisted  of  a  sub-set  of  these 
variables.  The  data  were  captured  to  a  text  file  that  was  then  opened  in  a  text  editor 
and  stripped  of  the  header  and  footer.  The  amended  text  file  was  then  loaded  into 
MATLAB®  software  for  detailed  analysis. 


Variable  ID  Col  Col#  Description  of  Variable 

in  Matrix 


1 

isid 

12 

1-12 

2 

mestim 

1 

13 

3 

idopp 

12 

14-25 

4 

icno 

12 

26-37 

5 

npclkb 

1 

38 

6 

rqxdkd 

1 

39 

7 

rqdat 

1 

40 

8 

rqilon 

1 

41 

9 

nkuah 

1 

42 

10:  npposn 

3 

43-45 

11;  iqivelc 

3 

46-48 

12:  gdop 

1 

49 

13:  sltdv 

12 

50-61 

14:  dtazi 

12 

62-73 

IS:  satlist 

5 

74-78 

16:  sltsta 

12 

79-90 

17:  ^isfix 

1 

91 

18:  iqll 

12 

92-103 

19;  npugsp 

1 

104 

20:  npuvsp 

1 

105 

21:  nknteh 

1 

106 

22;  nkntev 

1 

107 

23:  nkntep 

1 

108 

24:  nkprrs 

12 

109-120 

25:  nk^ 

12 

121-132 

26;  dfroln 

1 

133 

27:  cariock 

12 

134-145 

SV(PRN)  ID  for  raw  measurements  taken 

Local  time  of  week  (sec)  when  last  4  Hz  measurements  taken 

Measured  satellite  Doppler  frequency  (Hz) 

Carrier  to  noise  ratio  (C/No)  (dB-Hz) 

Receiver  clock  bias  (sec) 

Receiver  dock  drift  (seconds/sec) 

Current  position  latitude  (radians) 

Current  position  longitude  (radians) 

User  altitude  (m)  ^  \ 

Receiver  best  estimate  of  position  in  ECEF  coord  fiame  (m) 
Receiver  best  estimate  of  vdodty  in  ECtdl*  coord  fiame  (m/sec) 
Geometric  dilution  of  precision 
Satdlite  elevation 
Satdlite  aTimntli 

^^t<»iiitp<!  used  to  calculate  last  position  and  Vdodty  data 
Channel  status 

Fixing  status  of  receiver  (also  see  dfroln) 

Measurement  quality  of  pseudorange 
User  ground  speed  (m/s) 

User  vertical  speed  (m/s) 

F,sriinated  horizontal  accuracy  (m) 

Estimated  vertical  accuracy  (m) 

Estimated  position  accuracy  (m) 

TCalman  Filter  pseudorange  residuals  (m) 

Xalniaii  Filter  delta  range  residuals  (i^s) 

Flags  navigation  solution  as  differentially  derived 
Nuniber  of  measurements  without  cyde  slips  (4  Hz) 
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APPENDIX D-XR5-M  TECHNICAL  INFORMATION 


This  Appendix  describes  the  basic  technical  characteristics  and  configuration  of 
the  two  XR5-M  receivers  used  in  support  of  this  thesis.  Further  details  are  provided  in 
the  XR5-M  Control  and  Display  Unit  Software  Manual,  and  the  XR5-M  Data  Monitor 
Toolkit  Manual  [NAV93,  NAV96].  Navstar/Navsymm  Systems  Limited  is  a  subsidiary 
of  Telecom  Solutions  Europe  Ltd. 

The  XR5-M  is  a  twelve  channel  GPS  receiver  which  can  be  operated  in  either  a 
standalone  or  a  differential  mode.  The  unit  may  also  be  configured  as  a  base  station 
receiver  by  providing  real-time  differential  correction  data  to  mobile  GPS  receivers. 

The  receiver  may  also  be  configured  to  record  data  to  a  PC  for  subsequent  post¬ 
processing. 

The  Special  Committee  104  to  the  Radio  Technical  Commission  for  Maritime 
(RTCM)  Services  was  established  in  1983  to  develop  the  standards  for  differential  GPS 
corrections  [KAL86].  The  XR5-M  receivers  were  equipped  with  version  3.7  software 
which  operate  with  RTCM  Type  1,  Type  2,  Type  3  and  Type  9  Messages.  The  units  are 
specified  to  provide  differential  corrections  with  position  measurements  better  than  3 
meter  CEP.  The  latest  upgrade  to  the  software,  version  4.9,  supports  RTCM  message 
Types  20  and  21  and  provides  sub-meter  position  measurements. 

The  Control-Display  Unit  (CDU)  software  enables  an  operator  to  connect  the 
receiver  directly  to  a  PC  (using  MSDOS  version  3.3  or  higher)  in  order  to  initialize  and 
to  control  the  fimctions  of  the  GPS  receiver.  Communication  between  the  PC  and  the 
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XR5-M  is  through  either  the  COMl  or  COM2  port.  The  XR5-M  receiver  has  a 
conprehensive  list  of  variables  that  may  be  accessed  using  the  Data  Monitor  (DM) 
program.  A  master  symbol  file  of  all  variables  <ml237.sy>  must  first  be  loaded  into 
the  DM  program  before  a  subset  of  variables  can  be  generated.  A  detailed  description 
of  all  XR5-M  variables  is  provided  in  the  text  file  <datadic2.txt>  provided  with  the 
receiver  software.  The  master  symbol  file  must  also  be  loaded  each  time  the  DM 
program  is  used  in  order  to  execute  the  customized  symbol  variable  file.  The  master 
symbol  file  is  appropriately  named  for  the  number  of  receiver  channels  (12)  and  the 
software  version  (3.7).  The  DM  program  allows  an  operator  to  customize  the  data 
display,  and  to  generate  a  subset  of  variables  that  may  be  captured  to  a  data  file. 

Several  symbol  files  were  written  during  the  course  of  the  thesis  research.  The  list  of 
variables  used  for  the  collection  of  data  in  Test  12  is  provided  in  Appendix  C  and  the 
associated  symbol  file  name  was  <kdsymbB.sy>. 

During  the  preliminary  testing  phase,  the  first  attenpt  to  exit  fi-om  the  DM 
program  back  to  XR5-M  main  menu  using  the  command,  <$e>,  at  the  end  of  a  data 
capture  run  locked  up  the  PC  and  resulted  in  the  loss  of  the  data  file.  Subsequently 
Rising  the  command,  <$q>,  to  stop  data  collection  and  save  the  data  worked  fine  with  no 
lockups  over  the  six  weeks  of  testing.  This  procedure  was  limited  by  the  fact  that  it 
takes  the  user  back  to  DOS  pron:5)t  for  the  XR5-M  directory,  requiring  the  user  to  re-run 
the  XR5-M  executable  file  again  to  re-enter  into  the  XR5-M  main  user  menu. 

One  of  the  two  XR5-M  receivers  was  configured  as  a  base  station  and  was  used 
to  output  differential  corrections  to  the  second  receiver  configured  as  a  mobile  system. 
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Normally^  the  differential  corrections  from  the  XR5-M  base  station  receiver  would  be 
transmitted  via  a  Navsymm  DR5-96S  Radio  Modem,  The  RTCM  port  is  a  RS  232 
interface  connected  using  a  Series  723  12  Pin  Plug  as  shown  in  Figure  D.l. 


Figure  D.l  Series  723  12  Pin  Connector 

For  the  purpose  of  this  thesis  research,  the  RTCM  ports  of  each  of  the  receivers 
were  "hot-wired"  directly  together  to  enable  the  transfer  of  differential  corrections  from 
the  base  station  to  the  mobile  receiver.  The  base  station  receiver  RTCM  port  Pin  B 
(Radio  Control,  RC  also  known  as  Request  to  Send,  RTS)  was  wired  to  the  RTCM  port 
Pin  D  (Clear  to  Send,  CTS3)  on  the  mobile  receiver,  and  the  base  station  receiver 
RTCM  port  Pin  F  (Transmit  Data,  TXD3)  was  wired  to  the  RTCM  port  Pin  H  (Receive 
Data,  RXD3)  on  the  mobile  receiver  (see  Figure  D.l)  [NAV96]. 

Confirmation  of  differential  operation  of  the  mobile  receiver  was  verified  by 
several  methods.  The  receiver  reports  the  method  of  calculating  position  fixes  as  three 
dimensional  (3D),  two  dimensional  (2D),  and  no  valid  fix.  It  also  reports  whether  the 
position  fix  is  based  on  non-differential  or  differential  corrected  measurements.  The 
fixing  status  of  the  receiver  is  presented  on  the  Receiver  Measurement  Report  screen 
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(keystroke  "F3"),  ie.  three  dimensional  differential  (3D-DIF);  the  actual  differential 
corrections  sent  by  the  base  station  can  be  viewed  on  the  View  Differential  Corrections 
Report  screen  (keystroke  "F7")  and  the  differential  corrections  received  on  the  mobile 
receiver  can  be  viewed  using  the  same  keystroke.  Finally,  the  base  station  antenna 
position  (from  Appendix  B)  was  entered  as  a  waypoint  into  the  mobile  receiver  using 
the  Waypoint  Setup  Screen  (keystroke  "w").  The  difference  between  the  measured 
position  and  the  waypoint  measured  by  the  mobile  receiver  was  observed  on  the 
Waypoint  Navigation  screen  (keystroke  "n").  Prior  to  collection  of  data,  the  receiver 
was  verified  to  be  operating  with  a  valid  3D  differential  position  output. 

The  MSDOS  based  XR5-M  receiver  software  provides  access  to  a  very  large  set 
of  receiver  variables  which  can  be  downloaded  to  an  ASCII  file  and  then  opened  in  a 
text  based  word  processor  such  as  Microsoft^^  Word®  to  strip  off  header  and  footer 
lines  before  being  imported  into  MATLAB®.  All  automatic  features  such  as  spelling 
and  grammar  checking  should  be  disabled  to  ejqjedite  the  process  of  stripping  off  the 
header  and  footer  (given  that  some  of  the  file  sizes  were  in  excess  of 2000  text  pages/or 
15  Megabytes).  The  data  files  were  easily  read  into  MATLAB®  for  data  analysis. 
Overall,  the  numerous  variables  and  the  data  monitor  program  in  the  XR5-M  receiver 
make  it  an  excellent  receiver  for  education  purposes. 

The  XR5-M  volute  antenna  has  a  pre-att5)lifier  with  a  50  MHz  filter  centered  on 
LI.  Signals  outside  the  filter  bandwidth  are  reduced  by  40  dB  [BUT99].  The  pre-anp 
is  provided  due  to  the  cable  length  required  between  the  AFIT  building  rooftop  and  the 
navigation  laboratory.  The  volute  antenna  was  designed  for  the  dynamic  pitch  and  roll 
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of  a  maritime  environment.  For  future  testing  purposes  it  may  be  preferred  to  use  a 
choke  ring  antenna  in  order  to  reduce  multipath  effects.  A  single  antenna  shared  by  the 
two  receivers  may  also  reduce  the  multipath  effects  by  using  DGPS  since  the  multipath 
effects  would  be  common  to  both  receivers. 

The  XR5-M  has  two  intermediate  frequency  (IF)  stages;  the  first  IF  stage  is  at 
295  MHz  and  has  a  bandwidth  of  20  MHz,  the  second  IF  stage  is  at  24.58  MHz  and  has 
a  bandwidth  of  8  MHz.  The  receiver  uses  a  pre-correlation  A/D  and  it  has  a  2  bit 
quantizer  with  3  states  (above  a  positive  threshold,  below  a  negative  threshold,  and  in 
between  thresholds).  The  receiver  design  does  not  incorporate  a  frequency  lock  loop 
(FLL)  but  does  en:5)loy  a  phase  lock  loop  (PLL).  The  minimum  levels  of  C/No  (plus  or 
minus  a  few  dB-Hz)  for  receiver  operation  are  32  dB-Hz  for  acquisition  mode,  and  30 
dB-Hz  for  loss  of  lock  [BUT99]. 

The  variable  "nkntep"  provides  the  estimated  position  accuracy  (1  sigma)  and 
the  variable  "nknteh"  provides  a  measure  of  horizontal  position  accuracy.  A  navigation 
solution  is  considered  valid  if: 

i.  there  are  at  least  3  healthy  satellites  (for  a  2D  fix)  with  valid  ephemeris; 

ii.  the  estimated  accuracy  (1  sigma)  is  less  than  I  nautical  mile; 

iii.  the  velocity  is  less  than  515  m/s;  and, 

iv.  the  altitude  is  less  than  60,000  feet. 

For  the  installed  CDU  software  version  3.7,  GDOP  calculations  were  based  on 
the  five  satellites  used  in  the  navigation  solution.  The  latest  software  version  4.9  (not 
installed)  calculates  an  all-in-view  solution  and  the  DOPs  are  based  on  all  satellites  used 
in  the  solution. 
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Below  are  the  equipment  configuration  settings  recommended  for  anyone 
wishing  to  conduct  similar  testing  using  the  XR5-M  receivers.  From  the  DOS  prompt 
on  the  CDU,  change  directories  to:  >C:\cd  xr5ml2,  typing  <xr5m>  will  then  execute 
the  file  <xr5m.exe>  and  launch  the  main  CDU  program 

The  first  step  is  to  configure  the  two  receivers  via  the  CDU  of  each  PC.  Select 

keystroke  "1"  to  enter  the  Initialization  menu. 

The  following  inputs  were  used: 

Input  date  (d  m  y):  enter  current  date  ie.  17  12  1998 

UTC  time  (h  m  s):  enter  UTC  time(either  4  hrs  or  5  hrs  plus  local  time  for  AFIT) 

Initial  Lat:  enter  39  46  55. 128  N  and  it  will  automatically  display  it  as  39  46.9 1 88  N 
Initial  Long:  enter  084  04  59.648  W  and  it  automatically  display  it  as  084  04.9941  W 
Initial  Altitude:  enter  306.92  and  XR5-M  will  automatically  display  it  as  306.9 

The  operator  should  select  “Ctrl  s”  to  save  the  initialization  settings  to  file;  This 
research  effort  used  <kdinit>.  The  file  is  automatically  saved  in  the  DOS  directory 
>C:\xr5ml2\xr5m\data.  Any  future  initializations  can  be  ejq)edited  by  loading  this  file, 
simply  type  Ctrl  1  <filename>.  The  user  then  needs  only  to  amend  the  date  and  time  but 
initial  position  will  remain  valid  unless  antenna  location  is  changed  from  the  test 
location.  Once  the  parameters  are  correct,  the  operator  then  needs  to  select  “Ctrl  x”  to 
execute  the  settings.  The  XR5-M  will  reply  with  "Message  Acknowledged"  provided 
all  the  initialization  steps  are  con5)leted  correctly.  The  initialization  parameters  can  be 
checked  by  using  the  "FI "  function  key. 

The  Installation  Settings  are  slightly  different  for  the  two  receivers  depending  on 
whether  the  receiver  is  configured  as  a  base  station,  or  as  a  mobile  receiver.  Select 
keystroke  "2"  for  the  Installation  Screen. 
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Set  jR.eceivor  ^ 

Receiver  dynamics  set  to  5  for  DGPS  base  station  and  set  to  3  for  mobile  receiver 

Cable  Length  (m) 

Cable  Delay  (ns)  52.56 

Elevation  Mask  (deg)  ^ 

Timing  Pulse  Mode  Standard  1  PPS 

The  user  should  select  “Ctrl  s”  to  save  the  installation  settings  to  file.  This 


research  effort  used  <kdinstal>’  for  a  file  name.  This  file  is  automatically  saved  in 
>C:\xr5ml2\xr5m\data  directory.  The  user  needs  to  select  “Ctrl  x”  to  execute  the 


settings.  The  XR5-M  will  again  reply  with  "Message  Acknowledged". 

The  next  step  is  to  set  up  the  two  RTCM  ports  for  each  receiver.  The  RTCM 
settings  are  slightly  different  for  the  two  receivers  depending  on  whether  the  receiver  is 
configured  as  a  base  station,  or  as  a  mobile  receiver.  Select  keystroke  "5"  to  enter  the 


RTCM  Setup  Message  Screen.  Select  the  settings 


as  loiiows: 


Port  Function :  RTCM 

RTCM  Function;  select  "RTCM  Output  Single  Frame"  for  the  base  station  receiver 

RTCM  Function:  select  "RTCM  Input  &  Status  Output"  for  the  mobile  receiver 

Status  Message  Length:  Stan^d 

Lat/Lon 

Knots/Nm 
6  of  8 
4800 

n  •*  ODD 

Panty:  - 

Stop  Bits: 

Tx  Output  Interval: 

Tx  Output  Offset: 

Tx  Output  Offset: 

Authorization  Code 


Status  Message  Coordinates: 
Status  Message  Units: 
RTCM  format: 

Baud  Rate: 


1  s 
0  s 
100  ms 
0 


The  user  should  select  “Ctrl  s”  to  save  the  RTCM  settings  to  file.  This  research 


effort  used  <kddiffi>  for  a  file  name.  This  file  is  automaticaUy  saved  in  the 
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>C:\xr5ml2\xr5ni\data  directory.  The  user  needs  to  select  “Ctrl  x”  to  execute  the 
settings.  The  XR5-M  will  again  reply  with  "Message  Acknowledged". 

This  conpletes  the  initialization,  installation  and  RTCM  setup.  The  user  can 
then  proceed  to  the  Receiver  Measurement  Report  Screen  using  the  F3  function  key. 
Detailed  data  capture  is  conducted  using  the  Data  Monitor  program  Details  for  writing 
customized  symbol  variable  files  are  provided  in  the  CDU  Software  Manual  and  in  the 

Data  Monitor  Toolkit  Manual  [NAV93,  NAV96]. 

Overall  the  XR5-M  receivers  worked  very  well.  Graphical  presentation  of  real¬ 
time  data  using  the  XR5-M  software  was  limited  to  position  plots  only.  Ideally  real¬ 
time  graphical  plots  of  any  user-selected  parameter  would  be  desired,  ie.  Doppler 
tracking  frequencies  versus  time,  this  would  allow  for  the  assessment  of  jamming 
effectiveness  while  capturing  data.  Also,  saving  files  as  XR5-M  log  files  offer  a  very 
limited  ability  to  view  data  as  it  is  written  to  disk.  The  ability  for  the  user  to  pre¬ 
program  the  data  capture  period  in  the  Data  Monitor  program  prior  to  data  collection,  ie. 
two  hours,  would  be  a  desirable  feature  in  a  receiver  designed  for  testing  purposes. 


APPENDIX  E  -  HEWLETT  PACKARD  EQUIPMENT  DESCRIPTION 

The  three  main  Hewlett  Packard  (HP)  products  used  in  the  course  of  thesis  work 
were  an  HP  8643A  Signal  Generator,  an  HP  8563 A  Spectrum  Analyzer,  and  an  HP 
Variable  Attenuator  Model  #  3 94 A. 

The  HP  8643A  Signal  Generator,  Serial  Number-3221A00164  was  used  to 
generate  CW  and  swept  CW  signals.  The  signal  generator  operates  over  a  range  of 
0.252  to  2060  MHz  and  was  therefore  suitable  for  generating  jamming  signals  at  GPS 
frequencies.  The  output  was  set  at  -10  dBm  for  the  entire  test  period  and  the  HP 
Variable  Attenutator  was  used  to  adjust  jamming  levels.  The  accuracy  of  the  frequency 
output  by  the  signal  generator  was  verified  as  being  within  specification  by  the 
spectrum  analyzer;  however,  the  frequency  accuracy  of  the  signal  generator  was  0.375  x 
10"^  times  the  carrier  in  Hz  and  resulted  in  an  accuracy  of  ±590  Hz  [HEW92].  The 
XR5-M  receiver  was  used  to  measure  the  approximate  error  bias  in  the  signal  output  by 
the  fiignal  generator  by  observing  the  Doppler  frequency  offset  at  which  channels  were 
jammed.  This  bias  was  then  adjusted  for  using  the  frequency  selector  panel  on  the 
signal  generator.  This  was  discussed  fiirther  in  Section  4.9. 

The  HP  8563 A  option  026  Serial  Number  -  3222A02089  was  used  to  verify  that 
the  output  of  the  signal  generator  was  within  specifications.  The  HP  8563 A  Spectrum 
Analyzer  is  capable  of  measuring  signals  from  -1 19.9  dBm  to  +30  dBm  over  a 
frequency  range  of  9  kHz  to  26.5  GHz.  The  unit  was  warmed  up  for  at  least  five 
minutes  before  use.  The  accuracy  of  frequency  measurement  depends  on  several 
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factors  including  the  frequency  readout,  the  aging  of  the  analyzer,  temperature  stability, 
the  frequency  span,  the  resolution  bandwidth  plus  an  additional  10  Hz.  Based  on  the 
above  variables,  the  frequency  accuracy  of  the  spectrum  analyzer  was  estimated  at  ±460 
Hz  at  LI  [HEW91].  The  HP  8563A  provides  the  user  with  the  ability  to  program  and 
memorize  various  setups.  This  was  useful  for  setting  center  frequency,  frequency  span, 
and  resolution  bandwidth  to  best  measure  the  output  of  the  signal  generator.  The 
settings  could  be  recalled  at  a  later  date  for  various  tests. 

The  HP  Variable  Attenuator  Model  #  394A  is  adjusted  to  provide  6  to  120  dB  of 
attenuation  and  is  optimized  for  a  frequency  range  of  1000  to  2000  MHz.  Attenuator 
settings  of  120,  85,  80,  75,  70,  and  65  dB  were  used  dviring  the  tests.  Although  the  unit 
is  somewhat  dated  it  was  suitable  for  the  purpose  of  the  thesis  research. 

The  following  information  was  taken  from  technical  specification  sheet 
provided  by  HP  (the  sheet  was  retrieved  from  the  HP  archives  and  was  last  updated 
February  1, 1967  [HEW67]).  Absolute  accuracy  was  specified  as  ±1.25  dB  or  ±2.5 
percent  of  the  dial  reading,  whichever  was  greater.  This  would  accoimt  for 
approximately  ±2. 1  dB  of  error  for  an  attenuator  setting  of  85  dB.  The  attenuator 
handles  up  to  200  watts  average  depending  on  the  terminations.  Phase  shift  versus 
attenuation  is  negUgible  for  values  of  attenuation  above  approximately  20  to  30  dB 
because  the  imit  is  of  the  waveguide-beyond-cutoff  type  of  attenuator.  The  attenuator  is 
basically  free  of  phase  shift  with  attenuation.  Actual  phase  shift  results  from  changes  in 
input  and  output  reflection  coefficient  that  is  virtually  constant  for  attenuation  above 
approximately  20  to  30  dB.  During  data  capture  testing  the  attenuator  was  not  used 
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below  65  dB  attenuation.  The  reflection  coefficient  is  specified  as  less  than  0.23 
between  15  and  120  dB  of  attenuation  [HEW67]. 


The  unit  has  dimensions  of  140  mm  wide,  305  mm  long,  and  70  mm  deep  and 
weighs  6  lbs.  The  imit  had  a  list  price  in  1967  of  $550  US.  Hewlett  Packard  declared 
the  unit  obsolete  in  1984.  Newer  designed  models  have  subsequently  replaced  the 
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APPENDIX  F  -  MATLAB®  SOFTWARE  CODE 


Appendix  F  is  a  sample  listing  of  MATLAB®  software  code  written  for  three 
files  The  first  file  "appendt.m"  was  used  to  load  and  perform  calculations  on  the 
applicable  data  file  "dal7adec.txt"  for  Test  12  on  17  Dec  98  for  results  plotted  m 
Appendix  T.  The  associated  plotting  file  was  "mplott.m".  The  last  h^^AB  file 
"112ecef  m"  was  used  to  convert  fi-om  geodetic  latitude,  longitude  and  ellipsoid 

height  to  ECEF  coordinates. 


%  FUe  name  appendt.m 
%  Data  captured  on  17  December  98 
%  Times;  1430Zto  1705Z 
%  Sample  Rate:  One  sample  per  second 
%  Summary  of  Test 

%  This  was  a  data  run  with  10  periods  duplicated  for 
%  CW  at  70  dB  jamming  on  attenuator  (J/S  of  35.7)  for  ten  minutes  of 
%  jamming  followed  by  five  minute  periods  with  no  jamming 

.  . 

%  Note;  Correction  of  4.7  meters  in  ECEF  X  direction  applied 

%  Elevation  mask  angle  set  for  5  degrees  in  both  XR5  receivers 

• 

%  Load  captured  data  fi-om  text  file,  set  variables,  and  deteruMe  position 
%  errors.  A  list  describing  each  of  the  XR5  variables  is  provided  at 
%  Appendix  C.  The  names  of  XR5  variables  were  maintained  in  the 
%  following  Matlab  code, 
clear 

fid=fopen('dal  7adec.txt'); 

[data]=fecanf(fid.'%f, [145,9300]); 

fclose(fid); 

data=data'; 

isid=data(:,l;12); 

mestim=data(:,  13); 

idopp=data(;,  14:25); 

icno=data(:,26:37); 

npclkb=data(:  ,3  8); 

npclkd=data( :  ,3  9); 

nplat=data(:,40); 

nplon=data(:,4i); 

nkualt=data(;  ,42); 

npposn=data(;,43  ;45); 

npvelc=data(:, 46:48); 

gdop=data(;,49); 
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sltelv=data(:,50;61); 
sltazi=data(:, 62:73); 
satlist=data(:, 74:78); 
sltsta=data(:, 79:90); 
gpsfix=data(:,9 1); 
iqll=data(:,92:103); 
npugsp=data( : ,  1 04); 
npuvsp=data(:,  105); 
nknteh=data(:,  106); 
nkntev=data(:,  107); 
nkntep=data(:,  108); 
nlq)rrs=data(:,  109: 120); 
nkdprs=data(:,  121:132); 
dfsoln=data(:,  133); 
carlock=data(:,  134: 145); 

% 

format  long  e 

%  Load  "ALGN  MARK"  ECEF  coordinates  used  in  XR5  Base  Station  receiver 
afit=[5. 05979873 1 83428 le+05  -4.88229389975 19 16e-K)6 
4.059587157385 162e+06]; 

N=9300; 

fori=l:N 

%  calculate  norm  of  3D  miss  distance  w/o  correction 
miss=nonn((npposn(i,  1 :3)-afit)); 
missl=((*^PPOsn(i,  1 :3)-afit))'; 

%  write  3D  miss  distance  into  colunm  i 
misslb(:,i)=™issl ; 

%  write  ID  norm  of  miss  into  column  i 
missb(:,i)=miss; 
end 

%  set  time  for  x  axis  to  minutes 

t=(l:N); 

t=t'; 

t=t*l/60; 

%  transpose  of  matrix  to  get  N  by  1  for  missb 
%  this  is  uncorrected  position  error  vector 
missb=missb'; 

%  transpose  of  matrix  to  get  N  by  3  for  misslb 
%  this  is  uncorrected  position  error  in  x,y,  z  direction 
misslb=misslb'; 

%  create  new  matrix  of  3D  miss  distance 
miss2b=misslb; 

% 

%  subtract  out  actual  distance  between  GPS  antennas 
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%  this  was  an  approximation  in  ECEF  X  direction 
%  in  fact  the  4.7  meters  could  be  broken  down  more 
%  precisely  into  a  portion  in  the  x,y,2  direction 
% 

miss2b(:,  l)=misslb(:,  l)-4.7; 
fori=l:N 

%  determine  norm  of  miss  distance  after  correction 
%  for  distance  between  GPS  antenna 
miss3b(;,i)=norm(miss2b(i,  1 :3)); 
end 

%  transpose  of  corrected  miss  distance 
miss3b=miss3b'; 

% 

%  miss3b  is  the  variable  used  to  plot  3D  position  error  m  figures 

%  miss2b  is  the  variable  used  to  plot  position  error  in  X,Y,Z  direction  in  figures 


% 

%  mean  of  uncorrected  position  error 
meanmiss  1  b=mean(missl  b) 

%  mean  of  corrected  position  error 
meanmiss2b=mean(miss2b) 

%  mean  of  3D  position  error 
meanmiss3b=mean(miss3b) 

%  determine  2D  ( X  Y  )  norm  of  miss  distance  after  correction 
%  this  was  interest  only  for  2D  errors  not  including  height 
for  i=l:N 

niiss5b(:,i)=norm(miss2b(i,  1  ;2)); 
end 

niiss5b=miss5b'; 

meanniiss5b=mean(miss5b) 


% 

%  open  associated  plotting  program  mplott.m  to  see  Figures  for  test  12 
%  end  of  file 

%%%%%%%%%%%%%%%%%%%yo%%%%%%%%%'‘/o%%yo%%%%%% 
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The  MATLAB®  file  "mplott.m"  was  used  for  plotting  the  figures  used  in  Appendix 
T.  I  have  included  additional  comments  on  the  plotting  of  figures  in  Appendix  H. 
The  MATLAB*  code  for  mplott.m  is  provided  below. 


%  Appendix  T  Plots  for  Test  12  Results 
%  for  17  Dec  98  1430Z  to  1705Z 
%  Use  Print  position  in  Matlab  of  [.8  1  7  9.4] 

%  the  print  position  can  be  modified  to  optimize  plots/page  layout 

% 

%  plot  of  PRN  #  used  in  Nav  solution 
% 

figure(3) 

h=axes('Position',[0  0  1  l],'VisibleVo£F); 
axes('Position',[.2 .15  .7 .75]) 
plot(t,satlist);axis([0  155  0  31]) 
zoom 

ylabelCPRN  #  Used  In  Nav  Solution') 
xlabel('Time(minutes)') 
titleCStart  Time  1430Z  17  Dec  98') 
str(l)={' 

str(2)={'Figure  T.3.  PRN  #  Used  In  Nav  Solution  vs  Time'}; 

str(3)={' 

str(4)={"}; 

str(5)={'  T-3'}; 

str(6)={' 

set(gcf^'currentaxes',h) 

text(0.320,0.013,str,'FontSize',12) 

hold  off 

% 

%  Plot  of  PRN  #  assigned  to  each  tracking  channel 

% 

figure(4) 

h=axes(Tosition',[0  0  1  1], 'Visible', 'off); 
axes('Position',[.2  .15  .7  .75]) 
plot(t,isid);axis([0  155  0  31]) 


zoom 

ylabelCPRN  #  Assigned  to  Each  Tracking  Channel ) 

xlabel('Time(minutes)') 

titleCStart  Time  1430Z  17  Dec  98') 


str(2)={Ti^e  T.4.  PRN  #  Assigned  to  Each  Tracking  Channel  vs  Time'}; 


str(3)={' '}; 


str(4)={"}; 

str(5)={' 


T-4'}; 
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str(6)={'’}; 
set(gcf,'currentaxes’,h) 
text(0.25, 0.013, str,TontSize',  12) 
hold  off 
% 

%  Plot  of  GDOP  vs  time 
% 

figure(5) 

subplot(2,l,l) 

plot(t(l  :3000,;),gdop(l  ;3000,:)) 
title('Start  Time  1430Z  17  Dec  98') 
ylabel('GDOP') 
xlabel('Time(minutes)') 
subplot(2,l,2) 

plot(t(3000:9300,:),gdop(3000;9300,;)) 

ylabel('GDOP') 

xlabel('Time(minutes)') 

h=axes(Position',[0  0  1  l],'Visible',’off); 

axes(h) 

str(l)={"}; 

str(2)={' 

str(3)={'Figure  T.5.  GDOP  vs  Time'}; 

str(4)={' 

str(5)={' '}; 

str(6)={'  T-5'}; 

str(7)={"}; 
set(gcf'cuirentaxes',h) 
text(0.38,0.013,str,TontSize',12) 
hold  off 


% 

%  Plot  of  the  idopp  variable  vs  time 
%  over  the  entire  period 
% 

figure(6) 

h=axes('Position',[0  0  1  l],'Visible','off); 
axes('Position',[.2  .15  .7 .75]) 
plot(t,idopp) 


zoom 

ylabel(Doppler  Frequency  0£fset(Hz)  by  PRN  #') 

xlabel('Time(minutes)') 

titleCStart  Time  1430Z  17  Dec  98') 


str(2)={Tigure  T.6.  Doppler  Frequency  vs  Time  }; 


str(3H”}; 
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str(4)={'-}; 

str(5)={'  T-6'}; 

str(6)={"}; 

set(gcf,  'cuiTentaxes',h) 

text(0.35,0.013,str,TontSize',12) 

hold  off 

% 

%  this  is  detailed  plot  of  idopp  variable 
% 

figure(7) 

h=axes(Position',[0  0  1  l],'VisibleVoff); 
axes(Position',[.2  .15  .7 .75]) 
plot(t(l  :3000,  :),idopp(l  :3000, :)) 
zoom 

ylabel(T)oppler  Frequency  0£fset(Hz)  by  PRN  #') 

xlabel('Time(mmutes)') 

titleCStart  Time  1430Z  17  Dec  98') 

str(l)={"};  . 

str(2)={'Figure  T.7.  Doppler  Frequency  vs  Time  }; 

str(3)={' 

str(4)={' 

str(5)={’  T-7'}; 

str(6)={’ 

set(gcfi'currentaxes',h) 
text(0.35,0.013,str,TontSize',  12) 
hold  off 
% 

%  this  is  a  detmled  plot  of  idopp  variable 
% 

figure(8) 

h=axes(Position',[0  0  1  1], 'Visible', 'off); 
axes(Position',[.2  .15  .7 .75]) 
plot(t(3000:6000,:),idopp(3000:6000,:)) 
zoom 

ylabelCDoppler  Frequency  Ofifset(Hz)  by  PRN  #') 

xlabel('Time(minutes)') 

titleCStart  Time  1430Z  17  Dec  98') 

str(l)={' 

str(2)={Tigure  T.8.  Doppler  Frequency  vs  Time  }; 
str(3)={' 

str(4H"}; 

str(5)={'  T-8'); 

str(6)={' 

set(gct'currentaxes',h) 
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text(0.3  5,0.013,str,TontSize',  12) 
hold  off 
% 

%  this  detailed  plot  of  idopp  variable 
% 

figure(9) 

h=axes(’Position',[0  0  1  1], 'Visible', 'off); 
axes('Position',[.2  .15  .7  .75]) 
plot(t(6000:9300,:),idopp(6000:9300,:)) 


zoom 

ylabelCDoppler  Frequency  Offset(Hz)  by  PRN  #') 

xlabel('Time(minutes)') 

titleCStart  Time  1430Z  17  Dec  98’) 


strOH"}; 

st^(2)={Tig^re  T.9.  Doppler  Frequency  vs  Time ); 

str(3)={' 

str(4)={"}; 

str(5)={'  T-9'}; 

str(6)={"}; 
set(gcfi'currentaxes',h) 
text(0.35,0.013,str,TontSize',  12) 


hold  off 


% 

%  plot  for  gpsfix  variable  vs  time 
%  provides  indication  of  receiver  fix 
%  status,  2  means  3D  FDi,  1  means  2D  FIX 
%  and  a  zero  means  NO  FDC 
%  this  variable  can  be  used  with  variable  'dfsoln' 

%  to  show  whether  fix  is  based  on  a  differentially 
%  determined  solution  or  without  differential 
%  this  results  if  one  of  the  channels  used  in  the 
%  nav  solution  is  not  input  fi’om  the  base  station 
%  this  plot  shows  only  whether  it  is  2D,  3D  or  NO  FtX 
% 

figure(lO) 

h=axes('Position',[0  0  1  l],'Visible','off); 
axes('Position',[.2  .15  .7 .75]) 
plot(t,gpsfix);axis([0  155  -.5  2.5]) 


zoom 

ylabelCReceiver  Fix  Status') 
xlabel('Time(minutes)') 
fitle('Start  Time  1430Z  17  Dec  98') 
str(l)={"}; 

str(2)={'Figure  T.10.  Receiver  Fix  Status  vs  Time'}; 


F-7 


str(3)={"}; 
str(4)={'  ■}; 

str(5)={'  T-10'}; 

str(6)={”}; 
set(gcf^'currentaxes',h) 
text(0.35, 0.013, str,TontSize', 12) 
hold  off 
% 

%  plot  for  x,y,z  in  ECEF  coords  miss  distance 
%  in  meters  and  norm  ofx,y,z  in  meters 
% 

figure(ll) 

subplot(4,l,l) 

plot(t,miss2b(:,l)) 

titleCStart  Time  1430Z  17  Dec  98’) 

ylabel('X  axis  Miss(m)') 

subplot(4,l,2) 

plot(t,miss2b(:,2)) 

ylabel('Y  axis  Miss(m)') 

subplot(4,l,3) 

plot(t,miss2b(;,3)) 

ylabel('Z  axis  Miss(m)') 

xlabel('Time(minutes)') 

subplot(4,l,4) 

plot(t,miss3b) 

ylabel('3DMiss(m)') 

xlabel('Time(minutes)') 

h=axes('Position',[0  0  1  l],'Visible','off); 

axes(h) 

str(l)={"}; 

str(2)={Tigure  T.ll.  X,Y,Z  and  3D  Error(m)  vs  Time'}; 

str(3)={' '}; 

str(4)={’ 

str(5)={'  T-ll'}; 

str(6)={’ 

set(gcf;'currentaxes',h) 
text(0.30,0.002,str,'FontSize',  12) 
hold  off 
% 

%  det^ed  plot  for  x,y,z  in  ECEF  coords 
%  miss  distance  in  meters  and  norm  of  x,y,z  in  meters 

% 

figure(12) 

subplot(4,l,l) 
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plot(t(l  :8040,  :),imss2b(l  :8040, 1)) 
title('Start  Time  1430Z  17  Dec  98') 
ylabelCXaxisMissCm)') 
subplot(4,l,2) 

plot(t(l  :8040,  :),imss2b(l  :8040,2)) 
ylabelCY  axis  Miss(m)') 
subplot(4,l,3) 

plot(t(l  :8040,:),miss2b(l  :8040,3)) 
ylabel('Z  axis  Mss(m)') 
xlabel('Time(ininutesy) 
subplot(4,l,4) 

plot(t(  1 ; 8040,  :),miss3b(l  :8040, :)) 

ylabel('3DMiss(m)’) 

xlabel('Time(minutes)') 

h=axes(Position’,[0  0  1  l],'Visible’,'ofF); 

axes(h) 

str(l)={’ 

str(2)={'Figure  T.12.  X,Y,Z  and  3D  Error(m)  vs  Time'}; 

str(3)={' 

str(4)={' '}; 

str(5)={'  T-12'}; 

str(6)={"}; 
set(gcf^  'currentaxes'jh) 
text(0.30,0.002,str,TontSize',  12) 
hold  off 
% 

%  detailed  plot  for  x,y,z  in  ECEF  coords 
%  miss  distance  in  meters  and  norm  of  x,y,z  in  meters 
% 

figure(13) 

subplot(4,l,l) 

plot(t(8040:8340,:),miss2b(8040:8340,l)) 
title('Start  Time  1430Z  17  Dec  98') 
ylabelCX  axis  Miss(m)') 
subplot(4,l,2) 

plot(t(8040:8340,;),miss2b(8040:8340,2)) 
ylabelCY  axis  Miss(m)') 
subplot(4,l,3) 

plot(t(8040:8340,:),miss2b(8040:8340,3)) 
ylabel('Z  axis  Miss(m)') 
xlabel('Time(minutes)') 
subplot(4,l,4) 

plot(t(8040: 8340,  :),miss3b(8040;8340, :)) 
ylabel('3D  Miss(m)') 
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xlabel('Time(minutes)') 
h=axes('Position',[0  0  1  l],'Visible',’ofi’); 
axes(h) 
str(l)={' 

str(2)={'Figure  T.12.  X,Y,Z  and  3D  Error<m)  vs  Time'}; 

str(3)={’ 

str(4)={'  •}; 

str(5)={'  T-13'}; 

str(6)={"}; 
set(gcf^'currentaxes',h) 
text(0.30, 0.002, str,TontSize',  12) 
hold  off 
% 

%  detailed  plot  for  x,y,z  in  ECEF  coords 
%  miss  distance  in  meters  and  norm  of  x,y,z  in  meters 

% 

figure(14) 

subplot(4,l,l) 

plot(t(8340:9300,:),miss2b(8340;9300,l)) 
title('Start  Time  1430Z  17  Dec  98') 
ylabelCXaxisMissCm)') 
subplot(4,l,2) 

plot(t(8340:9300,:),miss2b(8340:9300,2)) 
ylabelCY  axis  Miss(m)') 
subplot(4,l,3) 

plot(t(8340:9300,;),miss2b(8340:9300,3)) 
ylabel('Z  axis  Miss(m)') 
xlabel('Time(minutes)') 
subplot(4,l,4) 

plot(t(8340;9300,:),miss3b(8340;9300,:)) 

ylabel(’3DMiss(m)') 

xlabel('Time(minutes)') 

h=axes(Position’,[0  0  1  l].'Visible',’ofif); 

axes(h) 

str(l)={"}; 

str(2)={'Figure  T.14.  X,Y,Z  and  3D  Error(m)  vs  Time'}; 
str(3)={' '}; 
str(4)={' '}; 

str(5)={'  T-14}; 

str(6)={' '}; 
set(gcf;'currentaxes',h) 
text(0.30, 0.002, str.TontSize',  12) 
hold  off 
% 
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%  this  is  a  plot  of  Receiver  Clock 
%  Bias  and  Receiver  Clock  Drift 


% 


figure(15) 

subplot(2,l,l) 

plot(t,npclkb) 

ylabel(Receiver  Clock  Bias  (seconds)') 

xlabel('Time(minutes)') 

title('Start  Time  1430Z  17  Dec  98') 

subplot(2,l,2) 

plot(t,npclkd) 

ylabel(Receiver  Clock  Drift  (seconds/second)') 
xlabel('Time(minutes)') 
h=axes(’Position’,[0  0  1  l],'Visible','ofi’); 


axes(h) 

str(l)={"};  .  .  n 

str(2)={'Figure  T.15.  Receiver  Clock  Bias  and  Drift  vs  Time }; 


str(3)={"}; 

str(4)={' 

str(5)={'  T-15'}; 

str(6)={"}; 

set(gc^'currentaxes’,h) 
text(0.26,0.002,str,'FontSize',  12) 
hold  off 


% 

%  this  is  a  plot  of  C/No  vs  time 
%  for  two  channels  per  page 
%  plot  of  receiver  channel  1  and  2 
% 

figure(16) 

subplot(2,l,l) 

plot(t,icno(:,l));axis([0  160  0  60]) 
ylabel('C/No  (dB-Hz)  Channel  T) 
xlabel('Time(minutes)') 
title('Start  Time  1430Z  17  Dec  98') 
subplot(2,l,2) 

plot(t,icno(:,2));axis([0  160  0  60]) 

ylabel('C/No  (dB-Hz)  Channel  2') 

xlabel('Time(minutes)') 

h=axes('Position',[0  0  1  l],'Visible','o£P); 

axes(h) 

str(l)={' 

str(2)={' 

str(3)={'Figure  T.16.  C/No  vs  Time'}; 
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str(4)={' 

str(5)={' 

str(6)={'  T-16'}; 

str(7)={"}; 

set(gc^  'currentaxes',h) 
text(0.38, 0.013, str,TontSize',  12) 
hold  off 
% 

%  plot  of  C/No  for  receiver  channel  3  and  4 
% 

figure(17) 

subplot(2,l,l) 

plot(t,icno(;,3));axis([0  160  0  60]) 
ylabel('C/No  (dB-Hz)  Channel  3') 
xlabel('Time(niinutes)') 
title('Start  Time  1430Z  17  Dec  98’) 
subplot(2,l,2) 

plot(t,icno(:,4));axis([0  160  0  60]) 

ylabel('C/No  (dB-Hz)  Channel  4') 

xlabel('Time(minutes)') 

h=axes('Position',[0  0  1  l],'Visible','off); 

axes(h) 

str(l)={"}; 

str(2)={' 

str(3)={Tigure  T.17.  C/No  vs  Time'}; 

str(4H"}; 

str(5)={’ 

str(6)={'  T-17'}; 

str(7)={"}; 

set(gcf^'currentaxes',h) 
text(0.38,0.013,str,TontSize',  12) 
hold  off 
% 

%  plot  of  C/No  for  receiver  channel  5  and  6 

figure(18) 

subplot(2,l,l) 

plot(t,icno(:,5));axis([0  160  0  60]) 
ylabelCC/No  (dB-Hz)  Channel  5') 
xlabel('Time(minutes)') 
titleCStart  Time  1430Z  17  Dec  98') 
subplot(2,l,2) 

plot(t,icno(:,6));axis([0  160  0  60]) 
ylabel('C/No  (dB-Hz)  Channel  6') 
xlabel('Time(minutes)') 
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h=axes('Position',[0  0  1  l],'VisibleVo£f); 

axes(h) 

str(l)={' 

str(2)={'-}; 

str(3)={Tigure  T.18.  C/No  vs  Time'}; 

str(4)={"}; 

str(5)={"}; 

str(6)={'  T-18'}; 

str(7)={"}; 

set(gc^'currentaxes',h) 
text(0.38, 0.013, str.'FontSize’, 12) 
hold  off 
% 

%  plot  of  C/No  for  receiver  channel  7  and  8 

figure(19) 

subplot(2,l,l) 

plot(t,icno(;,7));axis([0  160  0  60]) 
ylabel('C/No  (dB-Hz)  Channel  7') 
xlabel('Time(minutes)') 
title('Start  Time  1430Z  17  Dec  98') 
subplot(2,l,2) 

plot(t,icno(:,8));axis([0  160  0  60]) 

ylabel('C/No  (dB-Hz)  Channel  8') 

xlabel('Time(minutes)') 

h=axes(Tosition',[0  0  1  l],'Visible','off); 

axes(h) 

str(l)={' 

str(2H'’};  .  „ 

str(3)={Tigure  T.19.  C/No  vs  Time'}; 
str(4)={' '}; 
str(5)={' '}; 

str(6)={'  T-19'}; 

str(7)={"}; 

set(gcf;'currentaxes',h) 
text(0.38,0.013,str,'FontSize',12) 
hold  off 
% 

%  plot  of  C/No  for  receiver  channel  9  and  10 

figure(20) 

subplot(2,l,l) 

plot(t,icno(:,9));axis([0  160  0  60]) 
ylabel('C/No  (dB-Hz)  Channel  9') 
xlabel('Time(minutes)') 
title('Start  Time  1430Z  17  Dec  98') 
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subplot(2,l,2) 

plot(t,icno(:,10));axis([0  160  0  60]) 

ylabel(’C/No  (dB-Hz)  Channel  10’) 

xlabel('Time(niinutes)') 

h=axes('Position',[0  0  1  1], 'Visible', ’off); 

axes(h) 

str(l)={"}; 

str(2)={"}; 

str(3)={'Figure  T.20.  C/No  vs  Time'}; 

str(4)={"}; 

str(5)={’ 

str(6)={'  T-20’}; 

str(7)={"}; 

set(gcf,'currentaxes',h) 
text(0.38,0.013,str,'FontSize',  12) 
hold  off 
% 

%  plot  of  C/No  for  receiver  channel  1 1  and  12 

figure(21) 

subplot(2,l,l) 

plot(t,icno(;,ll));axis([0  160  0  60]) 
ylabel('C/No  (dB-Hz)  Channel  11') 
xlabel('Tune(minutes)') 
titleCStart  Time  1430Z  17  Dec  98') 
subplot(2,l,2) 

plot(t,icno(;,12));axis([0  160  0  60]) 
ylabel('C/No  (dB-Hz)  Channel  12') 
xlabel('Time(minutes)') 
h=axes('Position',[0  0  1  l],'Visible','off); 
axes(h) 

str(lH’'}; 

str(2)={' 

str(3)={'Figure  T.21.  C/No  vs  Time'}; 

str(4)={"}; 

str(5)={' '}; 

str(6)={'  T-21'}; 

str(7)={"}; 

set(gcf^'ciirrentaxes',h) 
text(0.38,0.013,str,TontSize',12) 
hold  off 
% 

%  this  is  plot  of  J/S  vs  time  using  "stairs"  function 

% 

figure(22) 
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ti=[  0  5  15  20  30  35  45  50  60  65 

75  80  90  95  105  110  120  125  135  140  150 

155]; 

y=[0  35.7  0  35.7  0  35.7  0  35.7  0  35.7  0 

35.7  0  35.7  0  35.7  0  35.7  0  35.7  0  0]; 

h=axes('Position',[0  0  1  1], 'Visible', ’off); 
axes(Tosition',[.2  .15  .7  .75]) 
stairs(tl,y);axis([0  155  0  40]) 
ylabel('J/S') 
xlabel('Time(inmutes)') 
title('Start  Time  1430Z  17  Dec  98') 
str(l)={"}; 

str(2)={'Figure  T.22.  Jamming  Level  J/S  vs  Time'); 
str(3)={' '}; 
str(4)={' '}; 

str(5)={'  T-22'}; 

str(6)={' '}; 
set(gcf, 'curTentaxes',h) 
text(0.33,0.0 13,str,TontSize',  12) 
hold  off 
%  End  of  File 

%%%%%%%%%%%%%%%%%%%%yo%%%%%%%%%%%°/o%%%%%/o 
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The  MATLAB®  file  "112ecef.m"  was  used  to  convert  geodetic  latitude,  longitude 
and  height  to  ECEF  coordinates  using  equations  fi-om  page  27  of  [KAP96]^  The 
file  was  specifically  written  to  use  the  coordinates  for  the  Doppler  Station  32058  in 
WGS  84  provided  in  Appendix  A  as  follows: 

Latitude  N  39  46  55.424 

Longitude  W  084  04  57.918 

EUipsoid  Height  27 1 .77  meters 

A  detailed  comparison  of  the  results  of  using  MATLAB®  and  Kaplan's  equations 
are  discussed  further  in  Appendix  B.  The  software  code  is  as  follows. 

%  File  name  is  "112ecef.m" 

%  This  file  was  generated  to  compare  lat/lon/ht  conversion  to  ECEF 
%  using  transformation  fi'om  Kaplan  page  27 

%  for  a  set  of  lat/lon/ht  and  ECEF  coordinates  given  by  1988  survey  for 
%  AFIT  Building  640  rooftop 
% 

%  Convert  geodetic  lat  long  ht  to  ECEF  for  Doppler  Antenna  location 
% 

b=63  56752.3 142; 
a=6378137; 
asq=a^2; 
bsq=b^2; 

e=sqrt(l-(bsq/asq)); 

latd=39*pi/180; 

latm=(46/60)'^pi/180; 

lats=(55.424/3600)*pi/180; 

lat=latd+latm+lats 

lond=84*pyi80; 

lonm=(04/60)’'‘pi/180; 

lons=(57.918/3600)*pi/180; 

lon=-(lond+lonm+lons) 

h=271.77 

xl=(a’'‘cos(lon)); 

x2=l+(l-e^2)'^((tan(lat)r2); 

x3=sqrt(x2); 

x=(xl/x3)+(h*((cos(lon))*cos(lat))); 

yl=(a*sin(lon)); 

y2=l+(l-e^2)*((tan(lat)r2); 

y3=sqrt(y2); 

y=(yl/y3)+(h’^((sin(lon))’^cos(lat))); 

zl=a*(l-e^2)*sin(lat); 

z2=sqrt(l-(e^2)*(sin(lat)r2); 
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z=(zl/z2)+(h*sin(lat)); 
u=[x  y  z]' 

% 

%  ECEF  values  WGS  84  from  survey  sheet  1988  for  Doppler  station  32058 
% 

ul=[506020.04  -4882282.48  4059593 .06]’ 

% 

%  Difference  between  Matlab/BCaplan  equations  and  DMA  calculated  position 
delta=u-ul 
% 

%  Difference  between  XR5  calculated  and  1988  Surveyed  location  for  Doppler 
%  station 


APPENDIX  G  -  SOURCES  OF  INTERFERENCE  TO  GPS 
RECEIVERS  AND  MITIGATION  TECHNIQUES 


SOURCES  OF  INTERFERENCE 

A  summary  of  the  various  forms  of  interference  to  the  GPS  signals  is  provided 
below.  For  more  detailed  information  the  reader  is  referred  to  reference  [PAR96]. 

TROPOSPHERIC/IONOSPHERIC 

Tropospheric/Ionospheric  interference  can  cause  several  effects  on  the  reception 
of  GPS  signals  especially  at  low  elevation  angles.  The  main  sources  of  ionospheric 
errors  are:  1)  group  delay  of  the  signal  modulation,  2)  carrier  phase  advance,  3)  Doppler 
shift,  4)  Faraday  rotation  on  linearly  polarized  signals,  5)  refraction  of  the  radio  waves, 
6)  distortion  of  pulse  waveforms,  7)  signal  ampUtude  fading  or  amplitude  scintillation; 
and  8)  phase  scintillations.  Tropospheric  effects  can  result  in  attenuation,  scintillation, 
and  signal  delay  [PAR96]. 

MULTIPATH 

Multipath  results  when  a  signal  arrives  at  a  receiver  via  multiple  paths  due  to 
reflection  and  diffraction  of  the  signal.  Differential  techniques  are  normally  used  to 
or  minimize  errors  that  are  common  between  two  receivers;  unfortunately, 
multipath  effects  can  be  very  different  between  a  Base  Station  and  a  Mobile  receiver. 
For  C/A  code  the  maximum  multipath  error  is  approximately  15  meters  [PAR96]. 
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HABMONICS  OR  INTER-MODULATION  PRODUCTS 

Radio  Frequency  Interference  (RFI)  can  have  a  narrow  bandwidth  of  about  25 
kHz.  One  usually  expects  harmonic  signals  to  be  weak  and  filtered;  however,  it  is 
possible  that  harmonics  or  products  fi’om  airborne  and  ground  based  transmitters  may 
fall  within  the  LI  fi-equency  band.  As  an  example  transmitters  at  787.71  MHz,  525. 14 
MHz,  315.084  MHz,  225.06  MHz,  105.028  MHz,  have  harmonics  at  LI  (the  2"‘‘,  3"^, 

5*  7*,  and  15*  respectively).  The  13*  harmonic  of  lower  VHF  fi-equencies  of  121.150, 
121.175, 121.200  MHz  and  the  12*  harmonic  of  131.200,  131.250,  and  131.300  MHz 
also  appear  within  the  GPS  bandwidth  [PAR96]. 


PULSED  INTERFERENCE 

Inadequately  filtered  signals  in  nearby  frequency  bands  may  generate  pulsed 
interference,  such  as  from  radar  systems.  In  the  case  of  intentional  interference,  a 
pulsed  jammer’s  power  is  cycled  on  and  off  and  the  pulse  repetition  frequency  and  pulse 
duration  modified  to  in^rove  jamming  effects  [PAR96]. 

IN-BAND  INTERFERENCE 

As  highlighted  with  the  incident  at  the  Air  Force  Research  Lab  [GWM98],  in- 
band  interference  can  be  caused  by  accidentally  transmitting  in  the  GPS  band.  Anyone 
intentionaUy  atten^ting  to  jam  the  reception  of  GPS  signals  would  choose  to  operate 
within  the  bandwidth  of  GPS  receivers.  A  wide  band  jammer  would  need  to  cover  a 
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sufficiently  wide  portion  of  the  spectrum  between  the  first  nulls  which  is  typically  at 
least  ±  1.023  MHz  for  C/A  code  receivers  at  LI  and  ±  10.23  MHz  for  P(Y)  code 
receivers  at  LI  and  L2  [PAR96]. 

OUT-OF-BAND  INTERFERENCE 

Out  of  band  interference  is  normally  caused  from  nearby  transmitters  and  poor 
RF  filtering  in  the  GPS  receiver,  although  satellite  communication  links  that  bracket 
the  GPS  signals  have  potential  to  cause  interference  problems.  As  the  demand  for 
valuable  frequency  spectrum  allocation  continues  to  grow  there  will  be  increased  efforts 
by  other  organizations  to  encroach  into  the  allocated  GPS  frequency  spectrum. 

SPOOFING 

Intelligent  spoofing  attempts  to  fool  a  GPS  receiver  by  creating  a  false  satellite 
signal  that  is  locked  on  to  by  the  GPS  receiver  as  though  it  were  a  vaUd  satellite  signal. 
The  information  in  the  false  signal  is  slowly  made  erroneous  and  the  GPS  receiver 
calculations  grow  in  error.  Although  one  would  normally  associate  the  generation  of 
PKN  codes  with  spoofing,  this  thesis  demonstrated  that  simple  CW  signals  can  be  used 
to  spoof  a  receiver  to  more  than  12  kilometers  (see  Section  4.8.1). 
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MITIGATION  TECHNIQUES 


INTRODUCTION 

The  two  general  forms  of  interference  suppression  may  be  categorized  as  spatial 
and  signal  processing.  Since  GPS  signals  are  transmitted  at  LI  and  L2  frequencies,  the 
signals  are  masked  by  obstacles  and  do  not  diffract  around  them  in  the  same  manner  as 
lower  frequency  signals.  If  the  interfering  source  is  ground  based  or  at  low  elevation 
angles  it  must  have  line-of-sight  with  the  intended  victim.  Signal  processing  techniques 
are  used  to  filter  out  interfering  sources  while  spatial  techniques  use  direction  of  arrival 
of  the  jamming  signals  to  reduce  the  jamming  effectiveness  while  enhancing  the 
direction  of  arrival  of  the  satellite  signals. 

RFAND  IF  FILTERING 

Out-of-band  interference  can  cause  non-linear  or  saturation  effects  in  the  low 
noise  amplifier  as  well  as  inter-modulation  effects  in  down-converter  mixers.  High 
peak  power  pulses  can  damage  the  sensitive  RF  front  end.  The  application  of  filtering 
at  RF  and  IF  reduces  these  effects.  A  C/A  code  receiver  may  have  a  RF/IF  filter 
bandwidth  as  wide  as  20  MHz.  This  bandwidth  may  be  narrowed  to  1-2  MHz  in  order 
to  gain  greater  selectivity  against  out-of-band  interference  [PAR96].  In  the  case  of  the 
XR5-M  receiver  the  first  IF  stage  is  at  295  MHz  and  has  a  bandwidth  of  20  MHz,  the 
second  stage  IF  is  at  24.58  MHz  and  has  a  bandwidth  of  8  MHz  [BUT99].  Rejection  of 
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narrow  band  interference  at  IF  may  be  achieved  by  using  adaptive  frequency  nulling 
filters. 

ADAPTIVE  INTERFERENCE  CANCELLER  (AIC)  2100 

Mayflower  Communications  Company,  Inc.  is  now  marketing  a  unit  titled  the 
“AIC  2100  Adaptive  Interference  Canceller  Module.”  The  unit  is  designed  to  operate  at 
RF  and  plug  directly  in  between  the  GPS  antenna  and  receiver.  It  is  reported  to  provide 
an  improvement  of  up  to  50  dB  of  suppression  to  narrowband  interference  for  C/A  and 
P(Y)  code  receivers  [UPT98].  It  should  be  mentioned  that  the  device  is  not  designed  to 
reduce  wideband  noise  jamming. 

ANTENNA  AND  ADAPTIVE  ARRAYS 

Several  techniques  are  used  in  a  GPS  receiver  to  take  advantage  of  the  fact  that 
the  interfering  sources  are  usually  earth  based.  A  simple  antenna  system  may  have  a 
typical  gain  of  0  dB  to  +  6  dB  for  the  upper  hemispherical  portion  but  may  drop  ofFnear 
the  horizon  (less  than  5  degrees  elevation)  to  approximately  -5dB.  For  an  elevation 
angle  less  than  5  degrees  other  error  sources  such  as  ionospheric,  atmospheric, 
multipath  and  physical  obstructions  may  be  more  predominant  than  RF  interference 
(except  for  the  case  where  the  RF  interference  is  intentional  i.e.  jamming). 

Spatial  interference  suppression  is  achieved  by  using  a  multi-element  antenna 
array.  Such  an  antenna  increases  the  gain  in  the  direction  of  the  various  satellite 
positions;  the  satellite  position  is  available  from  the  almanac  or  ephemeris  data. 
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Adaptive  weighting  of  antenna  elements  can  be  used  to  null  out  sources  of 
interference;  the  limiting  factor  is  the  number  of  elements  in  the  array  as  well  as  how 
close  in  azimuth  the  satellite  is  to  the  actual  interfering  source.  The  Controlled 
Reception  Pattern  Antenna  (CRPA)  is  an  example  of  such  a  system. 

At  the  US  Navy  test  range  at  China  Lake,  modeling  has  demonstrated  that  a  16 
kW  GPS  jammer  can  be  used  in  a  directional  antenna  beam  to  reduce  interference  with 
other  civilian  aircraft  at  locations  outside  the  main  beam  of  the  jamming  signal 
[BOG97].  Directional  jamming  may  be  used  in  the  future  as  civilians  become  more 
concerned  about  the  jamming  of  GPS  signals.  In  a  hostile  environment  the  operational 
benefits  of  jamming  certain  corridors  would  allow  fiiendly  forces  to  continue  to  operate 
in  a  GPS  jamming  environment. 

INTERFERENCE  SUPPRESSION  UNIT  (ISU) 

Electro-Radiation  Inc.  manufactures  the  ISU  that  has  been  tested  at  Holloman, 
AFB  with  four  ground  based  jammers  and  one  airborne  jammer  [BRASNY98].  The 
unit  uses  a  uniquely  designed  antenna  along  with  the  ISU  to  provide  an  additional 
margin  of  20  dB  for  broadband  noise  and  35  dB  for  narrowband  noise.  Future  testing  of 
this  device  in  the  AFIT  lab  environment  may  be  challenging  given  the  company’s  use  of 
polarization  as  a  means  of  reducing  interference. 
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AIDING  OF  THE  RECEIVER  TRACKING  LOOPS 

Once  interference  is  detected  it  may  be  possible  to  reduce  the  effects  of  the 
interference  by  narrowing  the  bandwidth  of  the  code  tracking  loop  or  the  earner 
tracldng  loop.  One  means  is  to  use  inputs  from  an  inertial  reference  unit  to  aid  the 
tracking  loops.  The  use  of  integrated  GPS  and  Inertial  Reference  Systems  is  likely  to 
increase  in  future  operations  due  to  the  complimentary  capabilities  of  both  systems. 

ANTI-SPOOFING  (A/S) 

Spoofing  of  satellite  signals  is  considered  to  be  challenging  given  that  the 
spoofer  normally  needs  to  know  the  receivers  position  and  velocity  in  order  to  create  the 
false  signal  that  will  have  the  correct  Doppler  shift  and  transmission  delay  of  the  codes. 
The  encryption  of  the  P  code  into  Y  code  is  designed  to  eliminate  this  threat.  The  C/A 
code  receiver  is  more  susceptible  to  spoofing  given  the  lack  of  anti-spoofing  capability, 
in  the  future  as  ground  navigation  aids  may  be  replaced  by  a  global  air  traffic  control 
system  using  GPS,  this  could  become  of  greater  concern  with  terrorist  activity.  This 
thesis  has  demonstrated  that  even  CW  interference  can  spoof  a  DGPS  C/A  code 
receiver  to  a  relatively  large  error.  As  a  result,  improvements  in  commercial  receiver 
design  to  reduce  CW  interference  will  serve  as  a  technique  to  reduce  the  impact  of  this 
simple  form  of  spoofing. 
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MULTIPATH 


One  method  of  reducing  multipath  interference  is  the  use  of  RF  absorbing 
materials  around  the  receiver  antenna.  Additionally  since  the  P(Y)  code  is  modulated  at 
ten  times  the  rate  of  C/A  code  the  P(Y)  code  chips  are  one  tenth  the  length  and  the 
multipath  errors  are  also  one  tenth  that  of  C/A  code.  Since  multipath  normally  arrives 
from  angles  near  the  horizon  the  effects  can  be  minimized  by  using  an  antenna  with  low 
gain  at  the  horizon.  A  choke  ring  antenna  is  commonly  used  to  reduce  multipath 
effects. 
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APPENDIX  H- SUMMARY  OF  LESSONS  LEARNED 


OVERVIEW 

A  summary  is  provided  in  Chapter  4  of  the  results  of  Tests  1  to  12.  The  purpose 
of  this  Appendix  was  to  capture  lessons  learned  during  the  course  of  thesis  research  for 
use  by  students  interested  in  foUow-on  research.  The  Appendix  is  divided  into  two 
main  areas  of  interest;  equipment  related  and  MATLAB®  comments. 

EQUIPMENT  RELATED 

The  following  comments  are  provided  with  respect  to  the  equipment  used  during 
the  course  of  thesis  research; 

a.  The  XR5-M  receiver  is  an  excellent  tool  for  education  purposes.  Access  to  a 
wide  range  of  variables  allows  a  user  to  tailor  data  capture  for  post 
processing  analysis; 

b.  The  XR5-M  data  monitor  program  requires  a  master  symbol  file 
<ml237.sy>  to  be  loaded  prior  to  writing  or  loading  customized  symbol 
variable  files.  Customized  symbol  variable  files  are  written  into  unique 
XR5-M  code  and  the  files  are  not  easily  modified  once  saved.  If  the  user 
wishes  to  add  a  new  variable  to  the  symbol  file  after  it  has  been  saved  it  is 
necessary  to  re-wiite  the  entire  set  of  desired  variables  to  file, 

c.  An  upgrade  of  the  XR5-M  receiver  firmware  to  version  4.9  would  provide 
the  capability  to  transmit/receive  RTCM  messages  Type  20  and  21.  These 
messages  provide  fiir  the  real-time  broadcast  of  differential  earner  phase 
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corrections.  Mobile  position  measurements  could  be  reduced  to  sub-meter 
accuracy  with  this  upgrade; 

d.  Mr.  Simon  Butcher  who  is  the  field  service  and  technical  support  manager 
fi-omNavstar  Systems  Ltd  in  England,  the  original  manufacturer  of  the 
XR5-M  receiver,  was  very  responsive  to  technical  questions  and  provided 
answers  promptly  and  to  the  level  of  detail  requested.  Mr  Butcher  may  be 
contacted  at  sbutcher@telecom.com  or  toll  free  firom  the  USA  at  1-800-456- 
0910; 

e.  The  maYimnm  data  rate  of  communication  between  the  two  receivers  via  the 
RTCM  port  was  4800  Baud.  The  maximum  rate  of  data  transfer  via  the 
Control  Display  Unit  (CDU)  was  9600  Baud.  The  CDU  9600  Baud  rate 
limits  the  amount  of  data  which  can  be  captured  fi'om  the  XR5-M  receiver 
each  second.  Adjusting  the  sample  rate  to  slower  periods  did  not  alleviate 
the  problem.  This  requires  judicious  selection  of  symbol  variables  for  data 
capture  if  post  processing  of  a  large  number  of  variables  is  desired, 

f  The  XR5-M  and  the  HP  signal  generator  and  spectrum  analyzer  should  be 
warmed  up  for  at  least  thirty  minutes  prior  to  use  for  detailed  measurements; 

g.  "Hotwiring"  the  two  XR5-M  receivers  together  to  provide  differential 
corrections  worked  well;  and, 

h.  The  HP  spectrum  analyzer  has  a  maximum  DC  input  of  0  volts.  The 

XR5-M  receiver  antenna  has  a  DC  powered  pre-amplifier  that  is  powered 
fi-om  the  XR5-M  receiver  through  the  antenna  cable.  There  was  potential  for 
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the  DC  output  from  the  XR5-M  receiver  to  feed  back  through  the  jamming 
cables  and  damage  the  spectrum  analyzer.  The  HP  variable  attenuator  was 
verified  that  it  blocked  DC  signals  to  the  spectrum  analyzer;  however,  a 
separate  DC  block  is  recommended  at  the  spectrum  analyzer  as  an  additional 
precaution  to  prevent  equipment  damage. 

MATLA^  COMMENTS 

On  both  PC  and  Mactintosh®  platforms  it  is  necessary  to  set  the  file  path  to  the 
desired  "m"  files.  Copying  of  MATLAB®  figures  into  Microsoft™  Word®  or 
PowerPoint®  resulted  in  degraded  quaUty  of  the  figures.  Efforts  to  edit  the  figure  once 
converted  to  MS  Word®  also  caused  undesired  changes  to  the  MATLAB®  figure. 

Saving  of  MATLAB®  figures  into  encapsulated  post-script  files  resulted  in  the  loss  of 
some  of  the  text  strings  which  had  been  added  to  the  MATLAB  figures.  The  best 
quaUty  was  achieved  by  doing  all  plots  in  MATLAB®.  Editing  and  annotating  specific 
figures  was  completed  using  the  "gtext"  MATLAB®  command  and  saving  the  particular 
figure  as  a  separate  "m"  file.  By  saving  the  figure  as  an  "m"  file  it  could  be  modified  by 
opening  the  actual  "m"  file  and  editing  the  text  code.  A  copy  of  all  MATLAB®  files 
(version  5.2)  on  compact  disk  (CD)  accompanies  the  AFIT  copy  of  the  thesis  should  the 
reader  desire  to  review  specific  test  results  in  detail.  The  MATLAB®  files  on  CD  are 
compatible  with  both  PC  and  Macintosh  version  5.2  of  MATLAB®  For  further  details 
the  reader  should  refer  to  the  <readme.doc>  file  included  with  the  CD. 
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As  a  rookie  MATLAB®  user  at  the  beginning  of  my  thesis  research  I  found  the 
following  MATLAB®  software  commands  to  be  useful: 

-  sptool  is  a  graphical  user  interface  tool  found  in  the  Signal  Processing 
Toolbox  that  was  used  extensively  in  conducting  analysis  of  data.  Variables 
may  be  imported  from  the  MATLAB®  workspace  and  saved  into  a  MATLAB® 

"spt"  file  format; 

-  fid(fopen)  fscanf  was  used  for  loading  the  data  text  file.  An  example  of  using 

these  commands  is  provided  in  Appendix  F; 

-  legend  vfas  used  with  a  -1  at  end  of  the  string  to  place  the  legend  off  of  the 

figure  plot; 

-  view  was  used  for  switching  and  rotating  a  polar  plot  of  satellite  azimuth  and 
elevation  to  a  compass  rose  with  west  on  left,  and  north  up; 

-  zoom  was  used  with  the  mouse  to  expand  a  portion  of  a  figure.  This  was 

beneficial  for  deciding  what  detailed  areas  to  plot; 

-  stairs  was  used  to  plot  J/S  levels;  and, 

-  gtext  was  useful  for  annotating  text  to  specific  locations  in  figures  using 

mouse/roUer  ball  input  device.  This  tool  was  used  to  label  the  various 
channels,  PRN  numbers  and  frequencies  on  the  Appendbc  I  to  T  figures. 

MATLAB®  "m"  files  did  not  execute  if  the  file  name  began  with  a  number. 
MATLAB®  Signal  Processing  Toolbox  "spt"  files  opened  with  "sptool"  if  the  file  name 
commenced  with  a  number. 
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Finally  running  MATLAB®  version  5.2  on  a  200  MHz  Pentium®  U  PC  would 
often  lock  up  the  PC  after  printing  several  pages  of  figures.  This  would  require 
rebooting  Windows  95®  when  using  "m"  files  with  large  matrices  (9000  rows  x  145 
columns).  This  problem  was  less  frequent  when  a  400  MHz  Pentium  n  PC  was  used. 
No  detailed  troubleshooting  of  this  problem  was  conducted;  however,  it  was  beheved  to 
be  associated  with  either  a  MATLAB®  memory  allocation  or  a  computer/printer 
network  problem. 
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Figure  P.8.  Doppler  Frequency  vs  Time 
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Figure  Q.7.  Doppler  Frequency  vs  Time 
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Figure  Q.15.  C/No  vs  Time 


Q-15 


C/No  (dB-Hz)  Channel  6  C/No  (dB-Hz)  Channel  5 


F 


StartTime1500Z06Dec  98 


60 


Time(minutes) 


60 


50 


40 


30 


20 


10 


20 


40 


60 

Time{minutes) 


80 


100 


120 


Figure  Q.16.  C/No  vs  Time 
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Figure  R.3.  PRN  #  Used  In  Nav  Solution  vs  Time 
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Figure  R.6.  Doppler  Frequency  vs  Time 
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Figure  R.7.  Doppler  Frequency  vs  Time 
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Figure  R.15.  C/No  vs  Time 
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Figure  R.17.  C/No  vs  Time 
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Figure  R.18.  Jamming  Level  J/S  vs  Time 
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Figure  S.6.  Doppler  Frequency  vs  Time 
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Figure  S.9.  Doppler  Frequency  vs  Time 
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Figure  T.6.  Doppler  Frequency  vs  Time 
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Figure  T.7.  Doppler  Frequency  vs  Time 
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Figure  T.9.  Doppler  Frequency  vs  Time 
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Figure  T.10,  Receiver  Fix  Status  vs  Time 
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Figure  T.12.  X,Y,Z  and  3D  Error(m)  vs  Time 
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Figure  T.16.  C/No  vs  Time 


T-16 


C/No  (dB-Hz)  Channel  4  C/No  (dB-Hz)  Channel  3 


60 


Start  Time  1430Z 17  Dec  98 


50  h 


40  h 


30 


20 


10 


|H|« 


I 

40 


20 


60  80  100 
Time(minutes) 


120  140 


60 1 - 1 - 1 - r 


20 

10 


■  I  _ I _ I _ 1—1 - l—l - L-iJ- 

20  40  60  80  100  120  140 

Time{minutes) 


Figure  T.17.  C/No  vs  Time 
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